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We consider static, spherically symmetric equilibrium configurations consisting of interacting dark 
matter/dark energy and embedded in an external, homogeneous chameleon scalar field. With the 
coupling function, form of which is taken to meet cosmological observations, we estimate the effect 
of such a nonminimal coupling on the properties of dark matter compact configurations. We show 
, that the masses and sizes of the resulting chameleon dark matter stars are smaller than those of 

systems with no field present. 
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I. INTRODUCTION 



Xj^ • During the past one and a half decades, a number of conclusive evidences have appeared for that the present 
Universe is expanding with some acceleration 0, Q . Such an acceleration could not be explained only by ordinary 
matter from which visible stars and components of galaxies in the Universe are made of. Invisible form of matter called 



I ' dark energy (DE), which works as a repulsive force against attracting gravity and drives the accelerated expansion, 
has been invoked. The observational data indicate that more than 70% of the energy density in the present Universe 
-£j ' could be assigned to DE. 

Another important notion used to explain the evolution of the present Universe is the gravitationally attractive 
__r invisible substance called dark matter (DM). Its contribution to the total energy density is estimated to be of order 
23%. DM is clustered on scales of the order of galaxies and clusters of galaxies, and its existence is so far evident only 
via its gravitational interaction. 

Despite the fact that the true origin of neither DM nor DE is currently known, various ways have been suggested 
to model them. The simplest approach is the so-called ACDM model, where DE is described by Einstein's A-term, 
and DM is supposed to be a pressureless fluid (cold dark matter). However, the well-known cosmological constant 
problem related to this model suggests one to look for other ways to describe the accelerated expansion of the Universe. 
Perhaps one of the most promising approaches to address the origin of DE are theories that include various fundamental 
fields \A . Other possible ways discussed in the literature are modified (non-Einstein) gravity theories [1, |3| and 
£T) • models with extra space dimensions || . 
7— I | Along these lines, cosmological models with DM and DE interacting with each other not only gravitationally but 
also by a direct coupling fMl^ (for a review, see, e.g., 0, G3) have a number of interesting features, and allow one 
. ' to solve some known problems, for example, the so-called coincidence problem. 

On the other hand, the question arises: if the interacting DM and DE do exist, how such interaction will affect the 
properties of small-scale objects consisting of such substance? In the bulk of the literature, various localized compact 
objects consisting of dark energy [l5rl2^ | and of dark matter [25rl29| have been studied. In these models, DE and 



DM interact either gravitationally or by direct coupling. The direct coupling has been considered in Ref. [27j, where 
a static compact configuration consisting of DM coupled to a scalar field responsible for DE was studied. In that 
paper, the assumption was that the mass of dark matter particles depends on the scalar field, that leads to significant 
changes in properties of the resulting configurations, such as their masses and sizes. 

In the present paper, we deal with another possibility when the mass of a scalar field which is responsible for 
the acceleration depends on the surrounding matter content. As the matter content, both ordinary (baryon) matter 
and dark matter can serve. The use of such scalar field, dubbed a "chameleon" scalar field, forms the grounds of 
the alternative model of dynamical dark energy which is proposed in Refs. [30T - I321 ] . With the assumption that the 
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chameleon scalar field interacts directly with dark matter, one can construct a cosmological model of the present 
Universe which is consistent with all current observations [10j. 

We use the model of Ref. 10] to describe compact objects consisting of dark matter which interacts directly with 
dark energy. We will identify the effect of this kind of nonminimal coupling on masses and sizes of resulting compact 
objects. 

The paper is organized as follows. In Sec. [11] the general set of equations is derived for equilibrium configurations 
consisting of interacting DM and DE. Choosing an equation of state of dark matter in the form of an ideal completely 
degenerate Fermi gas, in Sec. Ill Bl we write down equations for this particular case, and solve them numerically in 
Sec. IIIII In conclusion, we summarize the results obtained. 



II. EQUATIONS FOR EQUILIBRIUM CONFIGURATIONS 

We study the effect of the presence of direct interaction between a cosmological chameleon scalar field and dark 
matter on compact configurations consisting of such dark matter. To evaluate the effect, we specify the interaction 
by using the model of cosmological evolution of the Universe of Ref. [T(j- Our basic set up is that the dark matter, 
being embedded in an external, homogeneous chameleon scalar field, feels its presence not only gravitationally but also 
through the nonminimal coupling. Obviously, the characteristics of such mixed configurations will then be determined 
by the properties both of the dark matter and of the chameleon scalar field. 



A. Lagrangian and general set of equations 

General form of the Lagrangian of the system under consideration can be written as follows: 

L = -J^R + jd^V - V(<p) + /foOLDM • (!) 

Here ip is the real scalar field with the potential V((p), A = ±1 corresponds to the usual or phantom scalar field, 
respectively, Ldm is the Lagrangian of dark matter. The case / = 1 corresponds to absence of a direct coupling between 
dark matter and scalar field. However, even in this case the two sources are still coupled via gravity. Cosmological 
equations of Ref. [l(| can be obtained from the Lagrangian ([I) by taking A = +1 and £dm = where e denotes 
the energy density of dark matter. 

The Lagrangian ([1]) contains two functions of the scalar field: the potential energy V(ip) and the nonminimal 
coupling function /(<£>). Their form is specified proceeding from some general field-theoretical considerations in such a 
way as to be compatible with current astronomical and cosmological observations, and at the same time they should 
not contradict laboratory experiments [Tol. l30j - [3^ | . Following Ref. [l(j, we take (in units c= H= 1) 

f(cp) = exp [/% - ^ )/M P1 ], V(cp) = (M P i/cp) a , (2) 

where Mp\ is the Planck mass, a, /3, and ifo are positive quantities, with f} ~ O(l), and the mass scale is tuned to 
M v ~ 10 -3 eV in order for acceleration to occur at the present epoch. In this model, ifo corresponds to the current 
value of the cosmological scalar field. Note that the potential V(ip) is an example of a tracker potential in quintessence 
models (33[. Moreover, the above functions ([2]) are similar to those used within the framework of chameleon cosmology 
[3bTj32l | where the direct coupling between a cosmological scalar field and ordinary (baryon) matter takes place. 

In general, the presence of interaction between DM and DE may have potentially observable cosmological impli- 
cations |7Hl3j. On the other hand, at astrophysical scales, due to gravitational instabilities, interacting DM and 
DE may form compact dense objects composed primarily of dark matter. Clearly, physical properties of such mixed 
configurations will depend on the form both of the coupling function and of the potential energy of the scalar field. 
Using the particular form of interaction described by the Lagrangian ([T]) with the scalar functions determined by ^ , 
we clarify the question of how the presence of the nonminimal coupling influences masses and sizes of the resulting 
compact gravitating configuration. 

For this purpose, we start from a configuration consisting only of fermionic dark matter, with no scalar field present. 
Such objects can support themselves against gravitational contraction by the degeneracy pressure of fcrmions obeying 
Pauli principle, in the same way as in the case of neutron stars and white dwarfs. After introducing the nonminimal 
coupling to this system, we will track changes in its properties as they depend on the parameters appearing in the 
functions 

We will use Einstein's gravitational equations that allows to take into account relativistic effects. As a matter 
source in these equations, we take the energy-momentum tensor which is obtained by varying the matter part of the 
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Lagrangian (JTJ) with respect to the metric: 

T? =f[(e + P)u t u k - 8*p] + A d lV d k V - 6 k 



-d^tp-Viip) 



(3) 



where e and p are the energy density and the pressure of a dark matter fluid, respectively, and u % is the four- velocity. 
To describe a compact spherically symmetric configuration, we take static metric of the form 



ds 2 = e"d(x ) 2 - e x dr 2 



r 2 dQ 2 7 



(4) 



where v and A are functions of the radial coordinate r, the time coordinate x° = ct, and dVt 2 is the metric on the 
unit two-sphere. Using this metric and the energy- momentum tensor ([3]), we obtain three equations (hereafter we use 
c = h = 1): First and second are the (q) and (}) components of the Einstein equations, 
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and the third one follows from the law of conservation of energy and momentum, T k k = 0. Taking the i 
of this equation gives 



dp 
dr 
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(5) 
(6) 

1 component 
(7) 



In the above equations, the prime denotes differentiation with respect to r. Note that Eq. (JJj) differs from usual 
equation of hydrostatic equilibrium by the presence of an extra term in the right-hand side which is associated with 
the coupling function /. The appearance of this term is related to the choice Ldm = — £■ Another possible choice, 
L m = p, used in Refs. [23j, [2J, [3J, l35l| im plies that the term containing / in the right-hand side of Eq. is absent 
(on this subject see the discussion in |35j|). 

Equations ©-© must be supplemented by an equation for the scalar field which follows from the Lagrangian (JTJ: 
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Thus, we have five unknown functions is, A, tp, s, and p. Keeping in mind that e and p are related by an equation 
of state, there remains only four unknown functions. To determine these functions, we have four equations: the two 
Einstein equations ((SJ and ((BJ), the scalar-field equation and the equation of hydrostatic equilibrium (the 
modified Tolman-Oppcnhcimcr-Volkoff equation). 



B. The case of fermionic dark matter 



In order to perform calculations, we should choose some equation of state for dark matter. To elucidate how the 
interaction between DM and DE affects the characteristics of compact configurations, we will restrict ourselves to 
some simplest form of equation of state. Namely, we assume that dark matter is an ideal completely degenerate 
Fermi gas at zero temperature. Its equation of state is obtained using usual expressions for the energy density and 
pressure JHHi]: 
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(9) 
(10) 



Here my is the fermion mass, kF is the Fermi momentum, z = kr/mf is the relativity parameter, e and p are the 
dimensionless energy density and pressure expressed in units of m,j . 

In two limiting cases, this equation of state can be represented in simple power-law forms: (i) in the nonrelativistic 
case, z <§C 1, we get the polytropic law, p cx e 5 / 3 ; (ii) in the ultrarelativistic case, z 1, we have p = e/3. 
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We introduce dimensionless variables for the current radius r and mass M(r) of the configuration under consider- 
ation [Hj]: 



M r 
M=—, r=- with Ml 



Rl 



pi 



(11) 



Here Ml and Rl are characteristic mass and radius of a fermionic configuration obtained by Landau (for details, see 
Ref. [Hj]). Also, introducing the dimensionless scalar field, (f> = ip/Mpi, and rewriting the metric function A in terms 
of the Schwarzschild mass parameter M as e~ A = 1 — 2GM(r)/r 1 we represent Eqs. ((5|), (O, and ((8]) in the following 
dimensionless form: 
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For convenience, we drop hereafter the tilde. In obtaining Eq. ([12")) we have used Eq. ©. The functions / and V 
from ©, which appear in these equations, take the following dimensionless form: 



/(<£) = exp 0(0 - O )] , V(0) = <5(l/0) a with J = (M v /m / / 



(15) 



Thus, to describe the static configuration under consideration we have obtained three equations ([12|) - (fT4]) with the 
scalar field functions (ITSl). 



III. NUMERICAL RESULTS 



In this section, we solve the system of equations (fT2|) - (IT4l numerically for given a, j3, and S subject to the following 
boundary conditions in the vicinity of the center of the configuration r = 0, 

£-£c+|r 2 , Af~M 3 r 3 , 0~0 c + ^r 2 , (16) 

where (j> c and e c denote the central values of <j> and the energy density e at r = 0. The coefficients 62, M 3 , and cf>2 
can be found from Eqs. (|T2l - p^|) . In solving these equations, we will use the dimensionless equation of state which is 
given parametrically by Eqs. (|9]) and (|T0|) in terms of the variable z. 

We start the numerical procedure near the origin r w and proceed to the point r = r&, where the pressure is zero. 
We refer to the obtained solutions as internal solutions. The mass contained inside the sphere of the radius ri, will be 
treated as the mass of the dark matter star. 

As we suppose in Sec. [TTl our spherically symmetric configuration is taken to be embedded in an external, homoge- 
neously distributed cosmological scalar field tpo. To provide the smoothness of solutions along the radius, we require 
the internal solutions to sew solutions obtained for the region r > r& characterized by a nonzero scalar field energy 
density. Thus, for r > we proceed numerical solutions of equations (|13[) - (|14D retaining only the gravitational and 
scalar fields while the dark matter fluid is taken to be zero. At large distances, the scalar field reaches its asymptotic 
value (fo taken from [Toj . which can be rewritten in our dimensionless variables as follows: 



(17) 



Here fi^g and denote the dimensionless densities of dark energy and dark matter, respectively, measured in 

units of the current critical density. The factor l/^/Sn appears here since we take Mp\ = 1/y/G instead of the reduced 
Planck mass, Afpi = l/v / 8^rG, used in Ref. [Hj. Also, as in [Icj, we take the usual (nonphantom) scalar field, i.e., we 
set A = +1. 
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FIG. 1: The mass-radius relation for configurations with equation of state given by ((9))- (|10[) and functions /, V from (115[l . The 
parameter /3 — 0, 0.5, 1, 2, 2.3, 3, 5, 7, from top to bottom. The curve labeled by /3 = corresponds to configurations without a 
scalar field. 



Formally, at the points where the scalar field becomes equal to (j>o there is a nonzero gradient of the field. However, 
since for the values of the parameters a, (3, and S being used here (see below) the field goes to <fo only at r ^> rb, the 
magnitude of the gradient part of the energy density Tq is always much smaller than the potential energy V(<f>). This 
allows to consider the field in this region, to a good approximation, as homogeneous. 

Using the above procedure, we define the mass of the configuration under consideration as the mass of all matter 
(dark matter plus dark energy) contained inside the sphere of the radius r = n,. In Refs. (23l . l24j configurations 
consisting of ordinary (baryon) matter nonminimally coupled to a cosmological chameleon scalar field were called 
"chameleon stars." By analogy, let us call the configurations being considered here as chameleon dark matter stars, 
and denote their mass as Mcdms = Af(r&). 

Another physically relevant parameter is given by the coordinate f associated with the proper radius of the star. 
This quantity is defined such as to be invariant with respect to spatial coordinate transformations preserving the 
spherical symmetry, and can be presented as follows: 

f = f e x ' 2 dr' = [ 
Jo Jo 

Then, in dimensional variables, the proper radius is i£ pr0 p = f&i?£. Since Rl is equal to half the Schwarzschild radius, 
then we require F& > 2 to avoid black hole configurations. 

We turn now to a consideration of solutions obtained according to the above procedure. To perform numerical 
calculations, it is necessary to choose the values of the parameters a, j3, and 6. Here we fix a — 0.2, which has been 
used in the paper [l(| in describing the evolution of the Universe. The value of the parameter 6 appearing in (|15p 
depends both on the mass scale parameter M v (in Ref. [l(| it was taken as M v ~ 10 -3 eV) and on the fermion mass 
to/. Since at the moment it is not definitely known from which type of fermions gravitationally bound clumps of 
dark matter may consist of, various fermion particles are considered in the literature. This could be both su per l ight 
gravitinos with a mass of the order of 10~ 2 cV and superheavy WIMPs with a TeV mass scale [37j • Following [28j, we 
assume for definiteness that m/ lies in the range 1 eV < to/ < 10 2 GeV. Then the value of the parameter 5 is always 
considerably smaller than unity, and, as the numerical calculations indicate, the influence of the potential V(<f>) from 
Q15p on the solutions is negligibly small compared with other terms. I.e., the scalar field can be treated as having no 
potential energy on the scales associated with the fermion mass scale to/. 

Thus, we have only one free parameter (3 ~ O(l). Since its explicit value does not follow from any first princi- 
ples, we will vary its magnitude slightly, keeping track of changes in characteristics of compact configurations being 
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FIG. 2: Typical distributions of the dimensionless total energy density To from ((5|) (left panel) and of the current masses 
Mcdms = M(f) (right panel) as functions of the relative invariant radius f/fb for different /3. For both panels, the central 
dark matter energy density is taken as e c = 10~ 3 . In the right panel, the parameter ft runs the values 0, 1, 2, 3, 5, 7, from top 
to bottom. 



considered here. 

The results of numerical calculations of Eqs. (ir2|) - (Tl4l) with the boundary conditions (fT6)) are presented in Figs. Q] 
and[2] Figure [T] shows the mass-radius relation for different values of the parameter /3. In plotting the curves, we have 
used the values of the dimensionless central density of dark matter e c from (|16|) lying in the range 10~ 5 < e c < 10. 

In obtaining the solutions, we started from such central values <j) c at which the scalar field asymptotically, when 
r 3> ri,, went to the background value 4>$ from (| 1 T[) . The numerical calculations indicate that <f> c must be always less 
than that of given by Eq. (TP?)) . Since <po depends on j3, a situation may occur when for some (3 and s c it is necessary 
to set 4> c < 0. It is therefore inevitable in this case that </>(r) would cross zero at some point, giving a singularity in 
the potential V((f>), that is physically unacceptable. That is why, in plotting the curves in Fig. [TJ we restrict ourselves 
to only such values of e c which provide the positivity of (j>(r) everywhere along the radius. This implies that not 
all mass-radius curves have a maximum. There exists some critical value of j3 ~ 2.3 which still allows the curve to 
approach the maximum at some e c . At larger j3, the maximum is absent. 

In Fig. [TJ the curve labeled j3 = corresponds to the case studied in [28[ with no scalar field present. The points 
of the curve to the right of the maximum correspond to stable configurations: with increase of the mass of a star its 
radius decreases. The points of the curve to the left of the maximum refer to unstable configurations. The curves for 
configurations with the scalar field demonstrate the similar behavior for systems with j3 < 2.3 when the maximum is 
still present. At bigger values of j3, the absence of maxima means that all configurations are stable. 

By comparing the stable configurations with and without the scalar field, one can identify the changes brought 
about by the presence of the nonminimal coupling. Namely, with the increase of (3 that corresponds to a stronger 
interaction between dark matter and dark energy both the masses (at a fixed star radius) and the sizes (at a fixed 
star mass) of the systems with the field become smaller, as compared to the configurations without a scalar field. 
Moreover, the differences between masses and sizes in these two cases can be as high as the order of magnitude, at 
large values of (3 (cf . Fig. [5J) . 

These results agree with those we get when considering the radial distribution of the matter. The results of our 
calculations with some fixed central value of the dark matter energy density e c are shown in Fig. [2l As one can see from 
the plots presented in the left panel, configurations with the scalar field are characterized by a lower concentration 
of the matter towards the center. Also, taking into account the fact that the systems with the scalar field are more 
compact as compared to the field- free configurations having the same e c , we eventually obtain configurations whose 
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masses are less than that of the system without a scalar field, as demonstrated in the right panel of Fig. [2j 

Let us now derive an approximate formula for the dependence of the maximum mass of the stars on the parameter 
f3. To do this, first we convert the dimensionless variables (fTTj) to the physical ones. We remind also that not all 
mass-radius curves have a maximum. So, for the configurations with 5 > 2.3, we take the very left point of the 
mass-radius curve as corresponding to the maximum mass. Following [28[, we take the characteristic fermion mass 
scale asm/ = 1 GeV, for which the characteristic Landau mass Ml = 1.632A/©. Then the corresponding approximate 
expression for a maximum mass, M max = M max (J3), can be presented in the form 

M max « 0.627M S ( 1 ^) 2 e- - 1 ^ 3 . (19) 

This expression gives masses of the order of a stellar mass for the fermion mass to/ ~ 1 GeV. In the two extreme 
cases we have: (i) for superlight particles of mass m/ ~ 1 eV, we get very heavy and large configurations of mass 
M max ~ 10 17 — 10 18 M Q and the radius of the order of a typical galaxy cluster size, R ~ 10 24 cm; (ii) for superheavy 
fermions of mass to/ ~ 100 GeV, we get light and small objects of mass M max ~ 10~ 3 — 10~ 4 M Q and R ~ 10 2 cm. 

To summarize the results, we have studied the model describing compact gravitating configurations consisting of 
interacting dark matter and dark energy. Dark energy is modeled by the chameleon scalar field tp which has the 
nonminimal coupling to dark matter described by the function f(<p) in the Lagrangian ([T]). W e sp ecified the form 
of this function to provide agreement with current astronomical and cosmological observations [lCf . The particular 
case of a compact configuration with dark matter represented as an ideal completely degenerate Fermi gas has been 
studied in detail. In order to elucidate the role of the scalar field, we made a comparison of dark matter configurations 
supported only by the Fermi gas with the chameleon dark matter stars. The result is that the parameter j3 ~ 0(1) 
appearing in the coupling function /(</?) largely determines masses and sizes of the objects under consideration. In 
particular, as it is seen from Fig. [T] while j3 increases, both the masses (at a fixed star radius) and the sizes (at a fixed 
star mass) of the chameleon dark matter stars tend to become smaller than those of configurations without a scalar 
field. 
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